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We perform a systematic study to explore the accuracy with which the parameters of intermediate- 
mass black hole binary (IMBHB) systems can be measured from their gravitational wave (GW) 
signatures using second-generation GW detectors. We make use of the most recent reduced-order 
models containing inspiral, merger and ringdown signals of aligned-spin effective-one-body waveforms 
(SEOBNR) to significantly speed up the calculations. We explore the phenomenology of the 
measurement accuracies for binaries with total masses between 50 and 500 M@ and mass ratios 
between 0.1 and 1. We find that (i) at total masses below ~ 200 Mq, where the signal-to-noise-ratio 
is dominated by the inspiral portion of the signal, the chirp mass parameter can be accurately 
measured; (ii) at higher masses, the information content is dominated by the ringdown, and total 
mass is measured more accurately; (iii) the mass of the lower-mass companion is poorly estimated, 
especially at high total mass and more extreme mass ratios; (iv) spin cannot be accurately measured 
for our injection set with non-spinning components. Most importantly, we find that for binaries 
with non-spinning components at all values of the mass ratio in the considered range and at network 
signal-to-noise ratio of 15, analyzed with spin-aligned templates, the presence of an intermediate-mass 
black hole with mass > 100 Mg can be confirmed with 95% confidence in any binary that includes a 
component with a mass of 130 Mg or greater. 

PACS numbers: 04.30.Tv,97.60.Lf 

INTRODUCTION 

Advanced LIGO [T| and Virgo [2] detectors are expected 
to start taking data in late 2015 and 2016 [5], respec¬ 
tively. Compact binary coalescences are a key source of 
gravitational-wave (GW) signals for advanced detectors 
[e.g., 01E]- These may include binaries where one or both 
components are intermediate-mass black holes (IMBHs), 
with mass in the ~ 50 - few hundred Mq range. 

There is growing observational and theoretical evidence 
for the existence of IMBHs in globular clusters [see E, 
for review]. Observational evidence comes in the form 
of observations of ultra-luminous X-ray sources (ULX), 
cluster dynamics (though these are mostly sensitive to 
higher-mass IMBHs, whose GW signatures would be at 
frequencies below the detectors’ sensitive band), and, most 
recently, a tentative quasi-periodic oscillator observation 
of a 400 Mq IMBH [7j. On the theoretical side, a number 
of models have been predicted for IMBH growth, from 
direct collapse from very massive stars [e.g., [8] to runaway 
collision scenarios [a [ml or gradual growth through stellar- 
mass BH mergers [e.g., m or accretion US- 

Advanced GW detectors could observe inspirals of 
stellar-mass compact-objects into IMBHs in globular clus¬ 
ters [T3]. IMBH binary mergers are possible if the binary 
fraction in a globular cluster is sufficiently high to al¬ 
low the formation of two IMBHs m , or via mergers 
of two globular clusters with each other and subsequent 
coalescences of the IMBHs they host mm- Outside 
of globular clusters, merging compact binaries including 
IMBHs could form directly from isolated binaries com¬ 


posed of very massive (> 300 M 0 ) stars nzi- All of these 
scenarios could produce advanced-detector event rates of 
tens or more detections per year, though much lower rates 
are also possible. Meanwhile, future detectors with good 
low-frequency sensitivity, such as the proposed Einstein 
Telescope [T5], could observe up to thousands of IMBH 
binary mergers per year m and use high-redshift IMBH 
binary observations to probe the history of massive black 
hole formation mm- 

GW observations, which allow for a direct mass mea¬ 
surement, could provide the first completely unambiguous 
proof of the existence of IMBHs in the few-liundred-solar- 
mass range. If such IMBHs are discovered, their obser¬ 
vations would shed light on very massive star evolution 
and globular cluster dynamics. IMBHs could also prove 
to be particularly accurate probes of strong-field dynam¬ 
ical gravity, allowing for tests of the general theory of 
relativity [e.g., 1221 [23]. As the coalescence of IMBHBs 
is expected to be electromagnetically quiet, gravitational 
waves are likely to be the only means of observing these 
systems directly. For this reason, the LIGO and Virgo col¬ 
laborations have carried out searches for IMBH binaries 
in initial detector data (which did not yield detections) 
[241426] and intend to do the same in the advanced de¬ 
tector era with more sensitive instruments. However, in 
order to establish that an IMBH has been detected and to 
explore the scientific consequences of this discovery, it is 
necessary to analyze the GW signature of a coalescence in 
order to infer the parameters of the systems, particularly 
the component masses. In this paper we perform the 
first systematic study of the accuracy of IMBHB mass 
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measurements achievable with GW observations. 

Accurate models for GWs emitted from IMBH binaries 
must include the inspiral, merger, and ringdown phases of 
the coalescence. One of the most accurate available the¬ 
oretical waveform families are effective-one-body (EOB) 
models 07112171 . EOB is an analytical approach that 
combines post-Newtonian expansion, re-sunrmation tech¬ 
niques and perturbation theory with additional calibration 
of certain model parameters against waveforms obtained 
by numerically integrating Einstein’s equations on super¬ 
computers. These models are defined via a complicated 
set of ordinary differential equations in the time-domain 
and can be very computationally expensive to generate, 
limiting their use in parameter estimation studies to date. 
Novel reduced order modeling (ROM) techniques 0OH32] 
have allowed for the construction of fast and accurate sur¬ 
rogate models of EOB waveforms. In particular, we use 
the frequency domain ROMs |.‘i 11 [521 for EOB waveforms 
with spins aligned with the orbital angular momentum 
of the binary, defined in |33( 134] and implemented in 
LAL 1 35 1 . This allows us to perform simulations on a 
scale unprecedented for this class of sources. 



Figure 1: Characteristic amplitudes h c = \ff\h{f)\ of 
SEOBNRv2 05 injections with various total masses and 
mass ratios, at network SNR 15 used in this study. In 
black, the detector noise amplitude spectrum yjS n (f) of 
the Advanced LIGO design noise-curve (solid) and the 
Advanced VIRGO noise-curve (dashed). 


SIMULATIONS 

We performed a systematic study of the accuracy with 
which the masses and spins of the IMBHB could be re¬ 
covered from GW observations by Advanced LIGO and 
Advanced Virgo GW detectors operating at design sen¬ 
sitivity. For Advanced LIGO we used the zero-detuned, 
high-power detector configuration 05], and for Advanced 
Virgo a phenomenological fit to the design sensitivity 
curve EH, both of which are displayed in Figure [l] along 
with typical characteristic signal amplitudes. For our 
analysis we began generating the waveforms at a lower 
frequency of 10 Hz in both LIGO and Virgo instruments, 
allowing us to take full advantage of the low-frequency 
sensitivity of the instruments, which will be achieved to¬ 
ward the end of the decade. The use of zero-noise mock 
data sets to estimate parameter measurement accuracy 
relies on the assumption that the noise is stationary and 
Gaussian; although Berry et al. 05] demonstrated that 
non-stationary realistic noise does not significantly in¬ 
fluence parameter estimation for neutron star binaries, 
departures from stationarity (noise “glitches”) could play 
a larger role for low-frequency, short-duration IMBHB 
signals. 

To investigate parameter measurement accuracy, we 
analyzed a set of mock data sets (injections) with 
the LALInference [39] Bayesian parameter-estimation 
pipeline. This pipeline returns a set of samples from 
the joint posterior distribution for the signal parameters. 
We can readily convert this output into the innermost 
90% credible region, spanning from the 5th to the 95th 
percentile, on marginalized single-parameter posterior 


distributions, which we use as a proxy for measurement 
accuracy. We find this to be a more robust metric for 
measurement accuracy than the standard deviation of the 
highly non-Gaussian posteriors. 

We injected data with both SEOBNRvl and SEOB- 
NRv2 03] waveforms, and used corresponding single-spin 
ROM template families for recovery. We found result¬ 
ing measurement accuracies that are qualitatively and 
quantitatively similar, so we only show results from the 
more recent SEOBNRv2 model here. We included only 
the dominant l = m = 2 mode of the gravitational wave 
signal in our simulations, as this is the only mode included 
in the SEOBNRv2 reduced order model. This omission 
means that our results can only serve as a conservative 
estimate of the parameter estimation performance for 
IMBHB systems, since the higher frequency harmonics of 
the signal can carry information to further constrain the 
signal model, especially at high masses |40l - f42] . Further 
development of ROMs to include both higher harmonics 
and spin is necessary to provide timely results with the 
best possible accuracy. We did not include the cosmolog¬ 
ical redshift of the waveforms, so our results should be 
interpreted as measurements of the redshifted masses in 
the rest frame of the detectors. 

In order to explore measurement accuracy as a function 
of mass and mass ratio, we carried out injections for a 
broad range of total masses M to tai = m\+m 2 between 50 
and 500 Mq, where > m 2 are the component masses. 
For each mass, we injected systems at four mass ratios 
q = TO 2 /TO 1 of 1, 1/2, 1/4 and 1/10, always with non¬ 
spinning components. All simulated signals were oriented 









3 


such that the orbital angular momentum vector was in¬ 
clined at 30° to the line of sight between the Earth and 
the binary. Although the orientation and sky location 
were the same in all simulations, we do not expect the 
measurement of the mass and aligned spin parameters 
to be significantly affected by this choice, since we do 
not include higher modes which can couple mass ratio 
measurement to extrinsic parameter accuracy [42] . The 
distance was chosen to yield a constant coherent signal- 
to-noise ratio (SNR) of 15. 

For the Bayesian analysis, we used flat priors on 
the component masses within the range mi , m 2 € 
[5,1000] M 0 , limited the total mass to M tota i < 1000 M 0 , 
and limited the mass ratio m^/m-i > 0.01. We did not 
assume any of the source parameters were known when 
performing parameter estimation, allowing an isotropic 
prior on orientation, and a uniform-in-volume prior on 
binary location out to a luminosity distance of 15 Gpc (a 
redshift of ~ 1.9). The prior on the single aligned spin 
X was fixed to be flat in [—1,0.6], the range of validity 
of the SEOBNRvl [33] approximant. Since this prior 
distribution does not match the distribution of sources 
analyzed, we should anticipate that posteriors on indi¬ 
vidual injections can be centered away from the true val¬ 
ues, despite the self-consistency of LALInf erence, which 
has been demonstrated to produce X% credible intervals 
that contain the true value X% of the time [38] 139 ] [43 ] . 
For example, the low a a priori probability of high-mass 
extreme-mass ratio injections with non-spinning compo¬ 
nents, coupled to the asymmetry in the impact of remnant 
spin on the well-measured central frequency of the domi¬ 
nant ringdown harmonic [e.g., [44] , will lead to a typical 
over-estimate of the inferred total mass for such sources. 
This is compounded by the prior on distance p(d[/) oc d 2 L , 
which for a fixed amplitude tends to prefer higher mass 
sources at greater distances. 

RESULTS 

Mass measurement Figure [ 2 ] shows the 90% credible 
interval for the chirp mass, M. = rri\^ b rruJ 5 M t ot 1 , as a 
function of the total mass M to ta i- 

At lower masses, the signal is dominated by the ‘chirp¬ 
ing’ inspiral portion, and the phase evolution is a function 
of A4 at leading order, which is therefore the most strongly 
constrained parameter when analysing lower mass sys¬ 
tems g5j. We find that the width of the 90% credible 
interval on M is 0.3 — O.5M 0 at M to tai = 5OM 0 and 
0.7 — 3.5 Mq at M to t = 1OOM 0 . For comparison, the 
same interval is typically < 0.01 / 0.03 / O.1M 0 for binary 
neutron star systems / neutron star - black hole binaries / 
stellar-mass binary black holes, respectively [e.g., mats]- 

Meanwhile, as the mass increases, the inspiral moves 
to progressively lower frequencies and out of the sensitive 
band of the detector (see Figure [T]) and the merger and 



^total C^©) 


Figure 2: The 90% credible intervals for the chirp mass 
M. as a function of total mass M to t a i, for four mass 
ratios m 2 /mi- True values are indicated by the solid 
lines. As M tota i increases in the range 50 — 200 M 0 , the 
measurement of A4 becomes steadily worse as the 
inspiral potion of the signal moves out of the sensitive 
band of the detector. 


ringdown contribute an increasing fraction of the SNR 
(see Figure [3]). At masses above ~ 100 M 0 , the SNR 
is dominated by the merger and ringdown, and above 
~ 200 M 0 , by the ringdown. The ringdown frequency 
depends only on the total mass and spin of the merger 
product (the latter is a function of the mass ratio for 
non-spinning components). We therefore expect the total 
mass of high-mass systems to be better constrained than 
the chirp mass [this has previously been pointed out by 
Graff, Buonanno, and Sathyaprakash in ref. 148] : moreover, 
the accuracy of the M tota 1 measurement should become 
increasingly insensitive to the mass ratio. Indeed, this 
is the behavior we see in Fig. [4] which shows the 90% 
credible interval for the total mass. 

Alternatively, the mass measurement accuracy can be 
visualized by considering the 90% credible region widths 
on component masses. Figure [5] shows that component 
masses are generally harder to measure because of the 
significant uncertainty in the mass ratio typical for GW 
parameter estimation [45] . The mass of the lower-mass 
component m 2 is particularly poorly constrained, espe¬ 
cially at high masses and more extreme mass ratios, where 
only the total mass is encoded in the ringdown signature. 
The fractional uncertainty of the better-measured m\ 
component varies between 40% and 10% between M tota 1 
100 and 3OOM 0 . 

Spin All of our injections have non-spinning com¬ 
ponents, and the SEOBNRv2ROM waveform model 
which we used includes only a single spin parameter 
X = (m-ixi + m 2 X 2 ) /M, a combination of the dimen- 
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Figure 3: The relative SNR, the ratio of the SNR above 
and the SNR below the GW frequency at the innermost 
stable circular orbit (ISCO). We use the Schwarzschild 
ISCO /isco = 6 -3 / 2 /( 7 tM) which is strictly speaking 
only valid in the test particle limit. The relative SNR 
w.r.t. the ISCO depends on the mass-ratio. In contrast, 
the ratio of SNRs above and below half the ringdown 
frequency of each system is only weakly dependent on the 
mass ratio and reaches unity at a total mass of 150 Mg. 





Figure 4: The 90% credible intervals for M to tai- 


sionless spins \i = L ■ Si/mf aligned with the orbital 
angular momentum L that plays a dominant role in gov¬ 
erning the inspiral phase evolution through spin-orbit 
coupling ®[5S]- Figure [6] shows that the measurement 
accuracy of % decreases with total mass, as the inspiral 
moves out of the detector band. In general, % is not 
well constrained for non-spinning injections (x = 0), as \ 
values between ~ 0.2 and ~ — 0.5 are allowed, spanning 
about half of the prior range [—1, 0.6]. 


Measurability of parameters as a function of signal-to- 
noise ratio We also performed a series of simulations 
where we increased the signal-to-noise ratio from 5 to 100. 
The shape of the posterior probability density function 
approaches a multivariate Gaussian at high SNR; once this 
happens, we expect uncertainties on individual parameters 
to fall off as SNR -1 . As shown in Figure [ 7 ] this is indeed 
the case for SNRs larger than ~ 15. 


DISCUSSION 


A key question that will arise when a massive system is 
detected is whether we can confidently establish that the 
system contains an IMBH. As the coalescence of IMBHB 
systems is likely to be electromagnetically quiet, GW 
observations will be essential to measuring the parameters 
of these systems. Our results indicate that advanced GW 
detectors, using models which include inspiral, merger and 
ringdown, will be able to constrain the masses of detected 
IMBHBs, at least under the assumption of aligned spins. 
Figure [8] shows the 5% lower bound on the mass of the 
more massive component as that parameter increases. 
This indicates that, at a network SNR of 15 or greater, the 
accuracy of inference will be sufficient to determine at 95% 
confidence that a system with non-spinning components 
does indeed contain an intermediate mass black hole with 
mass > 100 Mgwhen the mass of at least one component 
is ~ 130 M 0 or greater. 

The SEOBNR reduced order models [TQ |35j have 
proven to be instrumental for performing systematic pa¬ 
rameter estimation studies with SEOBNR waveforms and 
provide speedups of up to several orders of magnitude. 
At the high total masses and low sampling rate used in 
this study time-domain SEOBNR waveforms are compar¬ 
atively quick to generate, but the speedup gained from 
ROM is still very significant. A single waveform evalu¬ 
ation with the ROMs is roughly 700 times faster than 
the likelihood computed from the time-domain SEOB- 
NRv2 model. Due to overhead costs the overall runtime 
of the simulations is roughly a factor 50 cheaper than with 
time-domain SEOBNR waveforms, allowing parameter 
estimation on sub-day timescales. As the waveform mod¬ 
els improve further we will be able to analyse a broader 
range of physical effects, including higher harmonics, with 
the rapidity the ROMs provide. 
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Figure 5: The 90% credible intervals for the component masses m i (left, larger companion) and m 2 (right, smaller 

companion). 
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Figure 6: The 90% credible interval for the effective 
dimensionless spin \- 


Figure 7: The width of the 90% credible intervals in 
chirp mass A4, total mass M to t a i, symmetric mass-ratio 
r) and spin y as a function of network SNR for 
mass-ratio q = 0.25 and total mass 100 M 0 . 
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